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I.  REVIEW  OF  GENERAL  THEORY 


This  section  is  not  meant  to  be  a  capsule  history  of  progress  in  super¬ 
conductivity  but  is  rather  an  attempt  to  give  some  background  for  the  theory 

being  applied  here  to  superconductivity  in  palladium  hydrides 0 

2 

In  1957  Bardeen,  Cooper  and  Schrieffer  (BCS)  developed  a  microscopic 
theory  of  superconductivity.  According  to  BCS  the  transition  temperature  is 
given  by 

Tc  =  1.14  <  w>  exp [-  1/N(0)V]  (1) 

where  ^  w  ^  is  a  typical  phonon  energy  and  N(0)V  is  the  interaction  strength; 
N(0)  is  the  electron  density  of  states  at  the  Fermi  surface  and  V  is  the  pair¬ 
ing  potential  arising  from  the  electron-phonon  interaction. 

It  should  be  kept  in  mind  that  Eq.  (1)  is  a  weak-coupling  approximation 


to 


N(0)V  / 


<  w>  /2kTr 


dx 


tanh 


(2) 


with  x  =  e/2kTc.  (Eq.  (2)  is  Eq.  (3,28)  of  Reference  2„) 
behavior  if  Eq.  (2)  is 


T 

c 


N(0)  V 
2 


(.  w  y 


The  asymptotic 


(3) 


Actually,  since  BCS  theory  is  a  weak-coupling  theory,  this  result  should  not 

3 

be  taken  too  seriously. 

Since  the  BCS  paper,  much  progress  has  been  made  in  understanding  the 
role  of  the  electron-phonon  interaction  in  normal  and  superconducting  metals. 

4 

Migdal  showed  that,  in  normal  metals,  the  electron-phonon  interaction  could 

1/2 

be  treated  accurately  even  for  strong  coupling  to  order  (m/M) 


1 
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Eliashberg  and  Nambu  extended  the  Migdal  treatment  to  the  superconducting 
state.  Eliashberg  Theory  takes  into  account  the  retarded  nature  of  the  phonon- 
induced  interaction  and  treats  properly  the  damping  of  the  excitations. 


n 

Based  on  the  Eliashberg  equations,  McMillan  made  an  extensive  study  of 
the  relation  between  microscopic  theory  and  observed  superconducting  transi¬ 
tion  temperature,.  The  centra^  result  of  McMillan ?s  paper  is  his  solution  of 
the  f inite-temperature  Eliashberg  theory  to  find  Tc  for  various  cases,  and 
the  construction  from  this  of  an  approximate  equation  relating  Tc  to  a  small 
number  of  simple  parameters.  (Evaluation  of  these  parameters,  however,  is 
not  at  all  simple.)  McMillan’s  Tc  equation  is 


Tc  =  (eD/1.45)exp(-  ) 

c  D  FA-y* (1+0. 62X) J 


(4) 


where  0D  is  the  Debye  temperature,  y*  is  an  effective  Coulomb  repulsion 
reduced  from  the  instantaneous  repulsion  y  by  the  fact  that  Coulomb  coupling 
is  propagated  much  more  readiLy  than  phonon  coupling,  and 

^  _  2j°°  dwq2(w)F(V) 


w 


(5) 


with  a2F(w)  being  the  electron-phonon  spectral  function,  F(w)  is  the  phonon 
density  of  states  and  a2(w)  represents  the  electron-phonon  interaction,  X  as 
defined  by  Eq.  (6)  is  a  dimensionless  measure  of  the  strength  of  a2F. 

Various  modifications  of  the  McMillan  equation  are  frequently  used: 

g 

Dynes  :  The  premultiplier  (6D/1.45)  is  replaced  by  <  w  >  /1.20. 


<  wn> 


2  1 

—  /  dw  a2(w)F(w)wn~ 

0 


(6) 


2 


oo 


This  is  the  "nth"  moment  of  g(w)  =  (2/Aw)a2F.  (nth  moment  =  /  dw  g(w)wn.) 

...  3  0 

This  is  the  notation  of  Allen  and  Dynes  and  is  consistent  with  McMillan. 
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Allen  and  Dynes  :  The  premultiplier  is  replaced  by 


F1  F-w-, 

1  l  log 

112 


OO 

where  w.  =  lim  w  =  expO^-  /  —  a2F(w)£nw);  F.  is  a  function  of  X  and  y*; 

108  n+0  A  0  w  1 

F 2  is  a  function  of  A,  y*,  and 


w 


log" 

Roughly  speaking,  X  represents  attraction  and  y*  represents  repulsion; 

bigger  X  leads  to  higher  T_;  bigger  y*  leads  to  lower  T  . 

u  c 

McMillan  has  shown  that  X  as  defined  by  Eq0  (5)  can  be  written 


X  .  NCO)  <I2> 

M  ^w2  y 


(7) 


where  M  is  the  ion  mass  and  I  is  the  electronic  matrix  element  of  the  change 
in  the  crystal  potential  U  as  one  atom  is  moved 

I(PP')  =  fa £  Cepp,  *  VUWp,df  (8) 


is  the  average  over  the  Fermi  surface  of  the  square  of  I„  While 
Eq.  (7)  i£  rigorous,  the  relationship  [Eq.  (4)  and  modifications]  between 
Tc  and  X  as  defined  by  Eq,  (5)  -  or  Eq.  (7)  is  not  on  such  rigorous  footing. 
Eq o  (4)  and  its  modifications  do,  however,  have  considerable  theoretical 
justification  (as  good  approximations)  and  considerable  experimental  verifi¬ 
cation  (again  as  good  approximations) . 
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One  of  the  standard  ways  to  obtain  X  (by  calculation)  is  to  use  the 

9 

Gaspari-Gyorffy  (GG)  approximation.  This  is  an  approximation  for  the 
numerator  (usually  called  q)  of  X  as  given  in  Eq.  (7)  and  expresses  q  in 
terms  of  scattering  phase  shifts  and  the  decomposed  electronic  density  of 
states  (decomposed  by  atomic  sites  in  real  space  and  by  angular  momentum 
label  £).  The  central  result  of  GG  is  to  approximate  n  =  N(0)<'l2^>  by 

~  2m  eF  y  2(A-f-l)sin2(6^.t.1-6^)nJj/njj/.fl 

N(Ep)  £  n/Un^if1)  ('9) 

where  N(Ep)  =  N(0),  the  6^  are  the  scattering  phase  shifts,  n^  is  the  density 

of  states  of  angular  momentum  l  at  Ep,  and  the  n^1^  are  the  "free  scatterer'! 

density  of  states.  (In  atomic  units  ~  is  replaced  by  1.)  Note:  see  also 

•ft 

Evans,  Gaspari,  and  Gyorffy,  reference  10. 

II.  SALIENT  EXPERIMENTAL  POINTS  IN  PALLADIUM  HYDRIDES 

The  first  known  enhancement  of  T£  due  to  H  was  Th4H15  with  Tc  of  9°K 
discovered  by  Satterthwaite  and  Toepke.11  (Pure  thorium  has  T£  =  1.37°K.) 
Palladium  Hydrides: 

1.  Pure  Palladium  is  not  superconducting  down  to  0.002°K  (ref.  12  -  this 

is  ref.  2  of  B.  Stritzker,  Phys.  Rev.  Lett.  42_,  1769  (1979). 

13 

2.  In  1972,  Skoskiewicz  found  superconductivity  in  PdHg  gy  with 

Tc  ~  4°K.  Skoskiewicz  used  an  equilibrium  technique  putting  the  Pd  in  H2  gas 
at  pressures  on  the  order  of  10  kbar. 

3.  Also  in  1972,  Stritzker  and  Buckel,^  using  an  implantation  technique 
achieved  T's  of  ~  9°K  and  11°K  in  PdH  and  PdD  respectively.  While  x  was 

X  X 


4 


not  known  accurately,  it  was  estimated  to  be  about  1.0  in  each  case.  The 
inverse  isotope  effect  (Tc  for  PdD  higher  than  Tc  for  PdH)  should  be  noted. 

It  appears  that  Pd/H  ratios  >^0.8  are  necessary  to  obtain  superconductivity 
(in  undamaged  Pd) „ 

Subsequent  to  Stritzker  and  Buckel,  the  Baranowski-Skoskiewicz  group 

15  16 

were  also  able  to  achieve  Tc’s  of  this  order  in  PdHx  and  PdDx.  * 

1  7 

4.  In  1974,  Stritzker  demonstrated  further  increases  in  Tc  upon 
addition  of  noble  metals.  Pd-Au-H,  T___„  =  13„6°K;  Pd-Ag-H,  T-_„v  =  15.6°K: 
Pd-Cu-H,  Tcmax  =  17°K,  The  addition  of  Rh  decreases  Tc„  It  appears  to  be 
accepted  as  fact  that  there  is  a  maximum  Tc  with  respect  to  the  noble  metal 
concentration  for  each  noble  metal.  The  "optimum"  Cu  concentration  is  45%, 

Ag  z  30%,  Au  16%.  Whether  there  is  a  maximum  Tc  with  respect  to  H  concen¬ 
tration  is  more  controversial.  Stritzker  states  that  there  is  but  it  must 
be  kept  in  mind  that  Stritzker  does  not  have  good  control  of  his  H  concen¬ 
tration. 

1 8 

5.  In  1979,  Stritzker  demonstrated  that  pure  Pd,  ordinarily  not  a 
superconductor,  becomes  superconducting  by  means  of  irradiation  at  low  temper¬ 
atures  with  He+  ions  with  a  maximum  Tc  of  3.2°K. 

18  19 

6.  For  Pd-Ni-H  systems,  Stritzker  and  also  Baranowski  and  Skoskiewicz 

find  an  initial  decrease  of  Tc  as  Ni  concentration  increases  above  10%.  The 
two  groups  differ  for  Ni  concentrations  above  10%.  Stritzker  finds  very  low 
Tc;  Skoskiewicz  finds  a  second  maximum  in  Tc  at  about  25%  Ni. 
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7.  Susceptibility:  For  both  PdAg  and  PdHv  (or  PdDv)  X  decreases 

AAA 

(roughly  linearly)  with  increasing  x  (see  references  1  and  21) .  For  PdAlx 

susceptibility  decreases  very  rapidly  with  x  (see  reference  21). 

21 

Comment:  While,  as  Ruvalds  states,  this  rapid  decrease  of  X  with  x 

makes  PdAlx  a  promising  candidate  for  superconductivity,  this  factor  is  no 

guarantee.  From  Ruvalds1  Figure  1,  the  X  vs.  x  curves  for  PdAg  and  PdH 

x  x 

are  quite  similar;  PdHx  (for  high  enough  x)  is  superconducting,  PdAgx  is  not. 

8.  Density  of  states:  The  total  density  of  states  at  Ep  drops  by  a 
factor  of  about  4  to  6  on  going  from  pure  Pd  to  PdH.  This  is  indicated  both 
by  experiment  [specific  heat,  susceptibility  -  see  references  in  Reference  1 
(page  259  for  specific  heat)  and  in  Reference  1A  (page  107  for  specific  heat 
and  for  susceptibility]  and  by  calculation  (Reference  1A  -  Switendick)  and 
Reference  20  -  Papaconstantopoulos,  et  al„ 

9.  Pd-Al-H:  Pd-In-H 

For  the  H  concentration  giving  momentum  Tc  for  a  given  A1  (or  In) 
concentration,  both  systems  show  a  very  slight  increase  in  Tc  as  A1  (In)  is 
added  initially  and  then  a  drop  in  Tc  -  this  drop  is  precipitous  in  A1  and 
somewhat  slower  in  In.^1 

III.  BAND  STRUCTURE  EVIDENCE  FROM  CALCULATION 

There  are  a  large  number  of  Pd,  PdH,  Pd-Ag-H,  etc.  band  structure 
calculations  in  the  literature  (see  references).  Some  earlier  calculations 
by  Switendick  are  very  fruitful  reading  (see  separate  report  by  me  (4/11/78) 
on  his  1970  Solid  State  Commun.  paper'*'^). 
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19A 

Switendick  pointed  out  the  fact  that  for  PdH,  using  names  like 

"proton"  or  "onion"  model  is  much  too  simple  -  these  names  have  connotations 

19B 

which  just  don't  hold  up.  He  also  pointed  out  that  the  H  s-electrons  fill 
three  distinct  types  of  states  (1)  the  H-induced  bonding  states  -  well  below 
Epj  (2)  the  0.36  hole  in  the  d  bands  of  Pd  (near  Ep) ;  (3)  The  sp  bands  of  Pd  - 
well  above  Ep  (see  Reference  30,  p.  146). 

Papaconstantopoulos  et  al  have  published  a  series  of  papers  on  PdH(D)x^0 
and  Pd-Ag-H  and  Pd-Rh-H  alloys^  (also  see  Miller  and  Satterthwaite^®) .  Two 
separate  memos  have  been  written  on  these  papers  (1978  -  The  comments  on 
Pd-Ag-H  refer  to  an  earlier  version  of  Reference  22).  These  are  APW  calcula¬ 
tions.  The  authors  use  the  Gaspari-Gyorffy  approximation  for  q.  The  main 
thrust  of  these  papers  is  that  the  s-like  density  of  states  at  the  hydrogen 
site  is  crucially  important.  The  overall  band  structure  is  in  reasonable 
agreement  with  Switendick ' s .  The  calculated  total  density  of  states  at  Ep 
shows  the  large  drop  on  going  from  pure  Pd  to  PdH  as  is  observed  experimentally. 

IV.  VARIOUS  MODELS 

The  various  models  constructed  to  explain  superconductivity  in  Pd  hydrides 
fall  into  three  categories:  (1)  Phenomenological,  (2)  Phenomenological  but 
include  attempts  to  put  realistic  quantitative  values  into  a  McMillan-like 
Tc  expression  or  into  some  modification  thereof,  (3)  Attempts  to  actually 
calculate  some  of  the  parameters  that  determine  Tc  from  band  structure 
calculation. 
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A.  Phenomenological 

23 

1.  Bennemann  and  Garland 

(a)  Spin-quenching:  Pure  Pd  metal  has  unfilled  d  states.  The  spin 
fluctuations  ("paramagnons")  associated  with  these  unfilled  d  states  destroy 
superconductivity.  In  the  BG  view,  the  main  function  of  H(D)  in  PdH(D)  is  to 
provide  s  electrons  which  tend  to  fill  these  d  states,  thus  quenching  the 
paramagnons;  coupling  between  acoustic  phonons  and  d  electrons  is  then 
sufficient  for  superconductivity. 

Comment:  Everyone  agrees  that  quenching  of  spin  is  necessary;  the 
question  is  whether  this  alone  is  sufficient.  I  feel  there  is  good  reason 
to  drop  the  BG  model  on  both  experimental  and  theoretical  grounds. 
Experimental:  Alloying  Pd  with  Ag  also  decreases  the  unfilled  d  states  but 
doesn't  lead  to  superconductivity  (see  ref.  33).  Theoretical:  Ganguly's 
1973  paper  showed  that  quenching  of  spin  alone,  even  under  the  most  favor¬ 
able  coulomb  condition  (y*  =  0,  i.e.,  complete  quenching)  leads  to  a  (A*-y*) 
which  is  not  large  enough  to  account  for  the  observed  T  . 

(b)  BG  Explanation  of  Isotope  Effect:  BG  explains  the  inverse  isotope 
effect  as  a  consequence  of  the  smaller  lattice  for  PdD  than  for  PdH;  they 
take  the  view  that  a  smaller  lattice  constant  is  better  for  higher  T  . 

Comment:  This  is  contradicted  by  the  observed25-27  negative  dTc/dP. 

One  should  note  that  this  negative  dTc/ dP  does  not  bear  on  the  BG  suppression- 
of-spin  idea. 
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2.  Auluck  :  The  point  of  view  taken  is  that  Pd  is  merely  a  host  for 
"metallic"  H. 

Comment:  This  is  contradicted  by  tunneling  measurements.  The  sizeable 

contribution  to  a2F(w)  from  Pd  indicates  that  superconductivity  isn’t  all 

from  H.  Skoskiewicz  et  al^  felt  that  the  observed  negative  dTc/dP  contra- 

24 

dieted  Auluck !s  model;  Rowe  ,  however,  feels  that  the  tunneling  data 

represent  the  strongest  contradiction  to  the  model.  Skoskiewicz  et  al1^ 

also  felt  that  the  observed  inverse  isotope  effect  contradicted  the  metallic 

hydrogen  model  -  I  suppose  on  grounds  that  some  sort  of  interaction  between 

Pd  and  H(D)  is  necessary  to  produce  this  inverse  isotope  effect, 

30 

3.  Miller  and  Satterthwaite  :  The  zero-point  motion  of  H  in  PdH  is 
larger  than  that  of  D  in  PdD.  These  authors  postulate  that  this  leads  to 

1  op 

differences  in  electronic  properties.  Utilizing  Switendick's  picture  of 
where  the  H  s  electrons  go,  MS  argue  that  the  large  vibrations  in  the  PdH  case 
cause  more  overlap  -  thus,  more  H  s  electrons  go  into  the  PdH  bonding  states, 
leaving  fewer  for  superconductivity. 

Comment:  This  may  be  consistent  with  Ganguly's  model  (below).  However, 
this  MS  model  will  give  rip  >  r|^.  This  difference,  combined  with  a  difference 
in  M  <w2>  could  well  lead  to  a  too-large  inverse  isotope  effect.  Calcula¬ 
tions  of  electronic  differences  between  PdH  and  PdD  is  beyond  the  present 
state-of-the-art  of  band  structure  calculation. 

B.  Phenomenological  Combined  with  Tc  Equations 
31 

1.  Hertel  :  Optical  phonons  don't  play  a  role. 
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Comment:  This  idea  seems  pretty  well  contradicted  by  the  tunneling 

measurements  of  Dynes  and  Garno32  for  PdH  and  of  Eichler  et  al33  for  PdH 

(see  also  Section  6063  of  Ref.  1)«  Hertelfs  model  is  also  inconsistent  with 

20  35  36 

the  calculations  of  Papaconstantopoulos  et  al.  ’  *  Hertel  purports  to 

show  that  if  optic  phonons  influenced  Tc,  a  regular  isotope  effect  would 

result;  I  am  not  able  to  follow  his  argument. 

Hertel' s  model  is  contradicted  by  the  coherent  neutron  scattering 

measurements  of  Rowe  et  al.  For  PcIDq  ^3  these  authors  found  large  dispersion 

in  the  longitudinal  optic  modes  (which  are  dominated  by  deuterium  motion) . 

They  also  found  that  the  frecuencies  of  the  acoustic  modes  of  PdDn  are 

U  #  o  j 

considerable  changed  from  those  in  pure  Pd,  again  in  contradiction  to  Hertel. 

These  authors  feel  that  a  simple  screened-pseudopotential  calculation  (such 

as  Hertel* s)  is  too  oversimplified  to  explain  the  superconducting  properties; 

calculation  of  the  phonon  dispersion  relation  in  Pd  hydrides  must  incorporate 

the  complete  electronic  band  structure  results  -  their  results  are  in  accord 
19B 

with  Switendick’s  band  structure  calculation.  (Rowe  et  al  consider  Hertel* s 
model  to  be  rigid-band.) 

Finally,  we  note  that  Hertel  has  X^  (in  PdH)  s  0o2  Xpd  (in  PdH)  (see 
his  Eq.  (3));  Papaconstantopoulos  et  al  (Table  II  of  Ref.  20)  have  XH  (in  PdH) 

S  2.1  Xpd  (in  PdH) c 

2.  Ganguly2^:  (see  also  References  37,  38)  Optic  phonons  play  a  major 

role. 

The  central  point  here  is  the  addition  of  XQpt  (X  optical)  to  Xacc 
(X  acoustical) o  Ganguly  adds  X  ^  by  means  of  a  three-square-well  model 
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(a  modification  of  the  McMillan^  equation  similar  to  that  of  Allender  et  al^ 
for  excitons)  and  gets  values  of  Tc  for  PdH  and  PdD  close  to  the  observed 
values  -  including  the  opposite  isotope  effect.  The  opposite  isotope  effect 
is  explained  by  the  increased  anharmonicity  of  H. 

*7  O 

Strong  support  for  Ganguly's  model  comes  from  the  tunneling  measurements0 
on  Al-oxide-PdH(D) .  These  measurements  showed  pronounced  structure  in  the 

current  voltage  curves  at  the  appropriate  place  for  H(D)  optical  phonons. 

33 

Eichler  et  al  consider  that  high  energy  local  modes  of  D  in  PdD 
definitely  share  in  the  electron-phonon  interaction  leading  to  superconduc¬ 
tivity;  this  conclusion  is  based  on  these  modes  showing  up  in  the  derivatives 
of  the  superconducting  tunneling  characteristic  curves. 

Ganguly's  ideas  are  supported  in  a  series  of  papers  by 
Papaconstantopoulos  et  al.  *  * 

Comment:  Ganguly's  model  is  consistent  with  most  of  the  observed  facts. 
One  open  question  is  associated  with  the  existence  of  an  optimum  H/metal 
ratio  for  maximum  T  .  Ganguly's  explanation  for  such  an  optimum  in  Pd/noble 
metal/H,  in  terms  of  "preferential  siting"  -  i.e.,  the  H  atoms  preferring  to 
sit  in  Pd-rich  regions,  remains  to  be  tested.  Also,  the  Ganguly  model 
predicts  that  the  inverse  isotope  effect  should  be  accentuated  in  Pd/Cu/H 
(smaller  "cage"  size);  this  does  not  seem  to  be  the  case  (see  Reference  1  and 
Figure  3  of  Reference  39)  but  this  is  not  really  clear-cut  evidence  against 
Ganguly's  model. 

While  the  views  of  Stritzker's  group  (see  pp.  3  and  4  of  Reference  44) 
are  similar  to  Ganguly's  regarding  phonon  modes,  I  believe  there  are  some 
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differences.  On  page  246  of  Reference  17,  Stritzker  says,  "The  difference 

in  the  maximum  Tc  values  of  the  three  Pd-noble  metal-H  systems  can  be 

described  by  the  isotope  effect •  ••",  Ganguly  would  say  that  force  constants 

are  involved  (and  not  just  mass  as  mass  enters  the  BCS  equation) . 

C.  Actual  Band  Structure  Calculation  Tc  Parameters. 

20  22  ^6 

1.  Papaconstantopoulos  et  al  *  *  *  -  These  authors  do  APW 

calculations. 


(a)  PdDx  and  PdHx  (primarily  from  Reference  20)  -  The  band  structure 
calculations  were  performed  self-consistently  for  two  choices  of  the  exchange 
parameters  within  the  Xct  scheme.  Spin- independent  relativistic  corrections 
(i.e.,  the  "scalar"  correcticns;  mass-velocity  and  Darwin)  were  included 
explicitly.  The  rigid-band  model  is  used  for  x  f  1.0.  From  the  band  struc¬ 
ture  calculation  they  compute  the  density  of  states  (decomposed  by  sites  in 
real  space  and  by  angular  momentum  quantum  number,  S,) .  From  the  decomposed 
density  of  states  they  obtain  r]  (the  numerator  of  X)  using  the  Gaspari- 
Gyorffy  approximation.9  A  central  feature  of  their  calculation  is  the 
importance  of  the  s-like  density  of  states  at  the  H(D)  site.  As  contrasted 
to  Miller  and  Satterthwaite,  the  present  authors  take  the  purely  electronic 
properties  of  PdH  and  PdD  to  be  identical,  thus  equals  r)p„  a2F(w)  for 
PdH(D)  is  taken  as  a2pdFacc(w)  for  0  <  w  <  Wj  and  as  a2H(-D)F0ptOO  for 

w2  —  w  —  00 *  M^w2^  (the  denominator  of  X)  is  taken  from  experiment. 40  The 

41 

Bennemann- Garland  expression  for  y*  is  used: 

0.26n(Ep) 


y*  = 


1  +  n(Ep) 


(10) 
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Using  these  parameters,  the  linearized  Eliashberg  equations^  are  solved 

following  the  formulation  of  Leavens.  They  have  also  obtained  solutions 

3 

for  Tc  using  the  Allen  and  Dynes  equation.  Their  calculated  results  for  Tc 
are  about  1°K  higher  using  Allen  and  Dynes  than  those  using  the  Eliashberg 
equations.  These  Tc  results  are  in  good  agreement  with  measured  values  and 
are  consistent  with  Ganguly's  model. 

Comment:  It  is  important  to  note  that  the  ratio 

M<w2>  H  =  1.2  M  <w2>d  (11) 

is  taken  from  experiment.40  With  such  a  ratio  one  is  well  on  the  way  to 
obtaining  the  inverse  isotope  effect.  I  feel  that  the  major  "confirmation" 
of  Ganguly's  model  by  these  authors  is  that  they  not  only  get  (with  no 
adjustable  parameters)  an  inverse  isotope  effect  but  also  obtain  values  agree¬ 
ing  with  the  experimental  values  both  as  to  individual  (PdH  and  PdD)  values 
of  Tc  and  as  to  the  relative  difference  in  T  . 

In  Reference  36,  the  Coherent-Potential  Approximation  (CPA)  is  used  for 
PdHx.  The  electron-phonon  interaction  and  Tc  as  a  function  of  x  are  in  good 

agreement  with  experiment  and  with  the  rigid-band  calculation.20  The  CPA 
36 

results  for  the  electronic-specific-heat  coefficient  as  a  function  of  x  are 
in  excellent  agreement  with  experiment  while  the  results20  of  the  rigid-band 
model  are  not.  The  CPA  results  are  consistent  with  the  central  conclusion  of 
Reference  20;  namely,  that  the  increase  of  Tc  with  x  is  mainly  due  to  the 
increase  of  the  hydrogen-site  electron-phonon  interaction  which  in  turn 
is  due  to  the  increase  with  x  of  the  s-like  density  of  states  at  the  H  site. 
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(b)  Pd-Ag-H  (primarily  Reference  22)  -  The  band  structure  calculation 

for  Pd,  Ag  H  (and  Pd,  Rh  H  )  was  performed  using  the  APW  method  and 
±  y  y  x  1  -y  y  x 

procedures  similar  to  that  of  Reference  20.  The  present  calculations  are 
not  carried  to  self-consistency.  The  virtual  crystal  approximation  (for  y) 

and  the  rigid-band  model  (for  x)  are  used.  Equation  (11)  is  used.  TQ  is 

3 

calculated  from  the  Allen  and  Dynes'  equation.  The  authors  feel  that  the 
assumptions  made  are  well  justified  within  the  ranges  of  x  and  y  used.  Results 
for  Tc  are  in  good  agreement  with  experiment.  One  adjustable  parameter  y*, 
is  now  used.  The  importance  of  the  s-like  density  of  states  at  the  H  site 
is  again  stressed. 

Comment:  It  is  important  to  note  that  the  "calculated"  maximum  Tc  with 

respect  to  Ag  concentration  is  not  pure  calculation  but  uses  an  empirical 
fact  as  input  -  namely,  the  assumption  that  one  can't  get  H  above  a  certain 
concentration  and  that  this  "maximum  H  content"  decreases  as  the  Ag  concen¬ 
tration  increases.  This  supports  ideas  of  Baranowski/Skoskiewicz .  Without 
the  assumption,  the  values  of  T^  in  Reference  22  would  simply  increase  as  y 
increases  (no  maximum) . 

We  note  that  for  Pdp  yAgg  ng^  appears  to  be  increasing  through 
x  =  1.0.  This  does  not  lend  any  support  to  getting  a  maximum  Tc  with  respect 
to  H  concentration.  ntota]_  i-s  decreasing  in  this  range.  (ns^  and  n^  for 
PdHx  behave  qualitatively  similarly  to  the  Pdg  yAgp  jHx  case.) 

Summary  of  models: 

1.  Spin  quenching  is  necessary  but  not  sufficient. 

2.  Optic  phonons  play  a  major  role. 
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3.  The  increased  anharmonicity  of  H  (over  D)  can  give  an  inverse 
isotope  effect  through 

M^w2^  =  cM^w2^ 

x  H  '  D 

with  c  =  1.2  (and  not  lo0)  and  with  =  r^. 

4^  The  calculations  of  Papaconstantopoulos  et  al  are  consistent  with 
Ganguly* s  model 0 

5.  Ganguly's  model,  while  leaving  some  questions  open,  seems  as  good 
a  starting  point  as  any, 

6,  The  Miller-Satterthwaite  idea  of  anharmonicity  contributing  to  an 
electronic  difference  in  PdH(D)  (and  to  ri^  >  cannot  yet  be  ruled  out. 

V.  PROPOSED  DIRECTIONS 
A.  Experimental 

(1)  Measurement  of  T  of  PdH  and  PdD  with  good  control  on  x.  Just 

V  X  A 

getting  a  good  H  profile  would  be  an  important  step.  Also  try  to  settle  the 
question  of  the  existence  of  a  maximum  Tc  with  respect  to  H  concentration. 
(The  papers  of  Stritzker  et  al  indicate  that  such  a  maximum  exists;  the 
Baranowski/Skoskiewicz  group  do  not  see  any  evidence  of  such  a  maximum  -  see 
page  332  of  Reference  15,  for  example.  For  definite  statements  by  Stritzker 
see  page  2  of  Reference  14,  Reference  1,  etc.) 

(2)  As  in  (1)  above  but  for  Pd^_yMyHx  (and  Pd^_yMyHx)  where  M  is  Cu,  Ag 
or  Au.  (Again  Stritzker  et  al  believe  that  a  maximum  Tc  with  respect  to  x 
exists  -  see,  for  example  page  403  of  Reference  45  and  page  262  of  Reference 
17.)  The  existence  of  a  maximum  Tc  with  respect  to  y  seems  pretty  well 
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established.  (See  also  Section  114  of  this  report.) 

(3)  Possible  preferential  occupancy  of  Pd-rich  regions  when  H  is  added 

38 

to,  say,  Pd^  ^Ag  «  This  is  an  idea  of  Ganguly?s  to  explain  the  (assumed) 
maximum  in  Tc  with  respect  to  H  concentration  in  Pd-noble  metal-H  systems. 
This  question  is  associated  with  (2)  above  but  is  also  of  interest  in  its  own 
right.  (Note:  If  there  really  is  ru^  maximum  Tc  with  respect  to  H  concen¬ 
tration  in  these  systems,  this  question  of  preferential  siting  becomes  less 
important.)  This,  if  it  exists,  would  be  one  way  in  which  increased  H  could 
lower  Tc  and  would  provide  a  competing  mechanism  to: 

(a)  Increased  H  more  optic  phonon  interactions  ->  higher  Tc« 

(b)  Increased  H  ->*  larger  lattice  higher  Tc<> 

Testing  this  idea  will  require  sets  of  phonon  spectrum  measurements.  I 
suggest,  in  Pd^^Ag^H^,  say,  keeping  x  fixed  (at  about  0.6,  say)  and 
increasing  y  -  see  where  spectrum  starts  to  change*  Ganguly’s  idea  requires 
that  fA  H  be  greater  than  fpd_H. 

(4)  Measurement  of  Tc  in  Pd  Ni  H  with  good  control  on  x  (and  y)  one 

i -y  y  x 

wants  that  x  which  gives  maximum  Tc  for  a  given  y.  (See  section  116  of  this 
report.)  I  don’t  know  if  the  Ni-H  force  constant  is  known.  The  Baranowski 

and  Skoskiewicz  curves  look  too  smooth  to  indicate  a  phase  change.  Can  Ni 

be  converting  some  of  its  own  s  electrons  into  d  electrons?  (Auluck28 

suggests  such  an  idea  for  PdH  -  addition  of  H  inducing  Pd  to  convert  some  of 

its  s-electrons  into  d-electrons.) 

In  my  opinion,  while  this  (possible)  second  maximum  in  Pd-Ni=H  is  very 
intriguing,  this  is  not  the  key  to  superconductivity  in  the  Pd  hydrides. 
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(5)  Try  Pd-In-H.  This  would  check  Ruvald's  measurements  but  now  with 
good  control  on  the  H  concentration. 

1  O 

(6)  Possibly  look  at  radiation  damage.  Stritzker  has  done  this  for 
pure  Pd  (see  Section  I I 5  of  this  report).  It  would  be  interesting  to  start 
with  superconducting  PdH,  say,  and  then  apply  radiation  damage,  noting  the 
effect  on  Tc. 

B.  Theoretical/Experimental 

Pursue  Ganguly's  idea  of  adding  another  element  to  "utilize  the  portion 
of  the  phonon  spectrum  between  the  acoustical  Pd  mode  and  optical  H  mode". 

This  third  element  would  (according  to  Ganguly)  have  to  form  its  own  sub¬ 
lattice  -  it's  not  sufficient  to  just  go  onto  the  Pd  lattice  substitutionally. 

C.  Calculation 

(1)  Pursue  an  idea  of  Stritzker  and  Luo^  regarding  the  observed  maximum 
in  Tc  with  respect  to  noble  metal  concentration.  They  suggest  that  increasing 
the  noble  metal  concentration  increases  the  screening  of  the  H  electrons  - 
this  should  cause  a  decrease  in  the  electron-phonon  interaction  (decrease  in 
X)  and  also  a  decrease  in  the  coulomb  repulsion  (decrease  in  y) .  The  former 
effect  tends  to  decrease  Tc,  the  latter  to  raise  Tc.  The  interplay  between 
these  two  effects  could  cause  the  observed  maximum.  I  have  done  some  very 
preliminary  and  purely  numerical  computations  on  this.^ 
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(2)  Continuation  of  band  structure  calculation  -  The  main  computer 
program  for  band  structure  calculation  (using  the  MPW  method)  of  PdHj  q  is 
"BSNACL"  this  program  does  not  yet  have  provision  for  relativistic  corrections. 
It  does  have  the  perturbation  procedure  for  computing  hydrostatic  AE.  It 
does  not  include  self-consistency  and  is  not  set  up  to  do  PdHx  for  x  ^  1.0. 

To  a  large  extent,  good  band  structure  calculations  for  PdHx  have 
already  been  done  by  Papaconstantopoulos  et  al.  The  feature  we  have  that  he 
does  not  have  is  the  ability  to  compute  small  AE  shifts.  Those  authors  would 
circumvent  this  by  calculating  both  E°  (for  lattice  constant  ag)  and  E'  (for 
lattice  constant  a^) .  This  is  probably  adequate  for  entities  like  total 
(and  decomposed)  density  of  states  at  Ep  (the  important  band  structure 
parameter  for  superconductivity  considerations) . 

They  probably  cannot  get  good  AE  values  for  specific  E(k)  due  to  their 
procedure  having  no  "fixed  zero"  of  potential  (see  pages  94  and  95  of 
Reference  47) . 
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